home *** CD-ROM | disk | FTP | other *** search
/ Merciful 1 / Merciful - Disc 1.iso / software / r / rexx_plus_compiler / rexxpluscompiler2.dms / in.adf / RexxPlPlot / Examples / examples.zoo / Example08.plot < prev    next >
Encoding:
Text File  |  1990-03-17  |  1.3 KB  |  51 lines
  1. /* Example08.plot - Demonstration of 3-d plotting */
  2. /* Modified for RexxPlPlot by Glenn M. Lewis - 9/12/89 */
  3.  
  4. address 'PlPlot'
  5. say 'Please be patient...'
  6.  
  7. NPTS=41
  8. opt.1=1; opt.2=2; opt.3=3; opt.4=3
  9. alt.1=60.0; alt.2=20.0; alt.3=60.0; alt.4=60.0
  10. az.1=30.0; az.2=60.0; az.3=120.0; az.4=160.0
  11.  
  12. title.1="\frPLPLOT Example 8 - Alt=60, Az=30, Opt=1"
  13. title.2="\frPLPLOT Example 8 - Alt=20, Az=60, Opt=2"
  14. title.3="\frPLPLOT Example 8 - Alt=60, Az=120, Opt=3"
  15. title.4="\frPLPLOT Example 8 - Alt=60, Az=160, Opt=3"
  16.  
  17.    do i=1 to NPTS
  18.       x.i = (i-1-(NPTS/2))/(NPTS/2)
  19.       y.i = (i-1-(NPTS/2))/(NPTS/2)
  20.    end i
  21.  
  22.    do i=1 to NPTS
  23.       xx = x.i
  24.       do j=1 to NPTS
  25.          yy = y.j
  26.          r = sqrt(xx*xx + yy*yy)
  27.          z.i.j = exp(-r*r) * cos(2.0*3.141592654*r)
  28.       end j
  29.    end i
  30.  
  31.    'plstar(1,1);'
  32.  
  33.    do k=1 to 4
  34.       'pladv(0);'
  35.       'plvpor(0.0,1.0,0.0,0.9);'
  36.       'plwind(-1.0,1.0,-0.9,1.1);'
  37.  
  38.       'plw3d(1.0,1.0,1.0,-1.0,1.0,-1.0,1.0,-1.0,1.0,alt.k,az.k);'
  39.       'plbox3("bnstu","x axis",0.0,0,"bnstu","y axis",0.0,0,"bcmnstuv","z axis",0.0,0);'
  40.       'plot3d(x,y,z,work,NPTS,NPTS,NPTS,opt.k);'
  41.       if k=2 | k=4 then 'plside3(x,y,z,NPTS,NPTS,NPTS,opt.k);'
  42.       if k=3 | k=4 then 'plgrid3(0.0);'
  43.       'plmtex("t",1.0,0.5,0.5,title.k);'
  44.    end k
  45.  
  46.    'pltext();'
  47.    'plend();'
  48.  
  49. exit 0
  50.  
  51.